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Abstract 

We have considered the two-spin interaction spherical spin-glass model with asymmet- 
ric bonds (coupling constants). Besides the usual interactions between spins and bonds 
and between the spins and a thermostat with temperature T a there is also an additional 
factor: the bonds are not assumed random a priori but interact with some other thermo- 
stat at the temperature Tj. We show that when the bonds are frozen with respect to the 
spins a first order phase transition to a spin-glass phase occurs, and the temperature of 
this transition tends to zero if Tj is large. Our analytical results show that a spin-glass 
phase can exist in mean-field models with nonrelaxational dynamics. 
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Investigation of nonequilibrium stationary states for systems without detailed balance 
is an important problem. Systems of this type are widely used in nonequilibrium ther- 
modynamics 0, ||, biology and biophysics ||, theory of artificial neural networks f|, etc. 

If detailed balance (or potential conditions) holds, then the stationary state of a 
macroscopic system can be described by the Gibbs distribution (maybe with some gen- 
eralizations), which is independent from details of the dynamics @, 0. In the oppo- 
site case such a universal and simple distribution does not exist, and the situation is 
strongly dependent on the details of dynamics and the concrete form of detailed balance 
violation. In general, we have in a stationary state without detailed balance only time- 
translation invariance || [/]]: average multitime quantities C(ti,..,t n ) (time ordering is 
assumed t\ > ... > t n ) depend only on the corresponding time differences t k — t h k < I 
(for example, C{t) = const, C(ti,t 2 ) = C{t\ — t 2 )r--) 

In the present paper we consider the case when the detailed balance condition is 
violated by nonconservative forces of some specific type: the mean-field spherical spin- 
glass model with asymmetric bonds (coupling constants). This type of nonconservative 
forces was inspired by applications in the theory of artificial neural networks M , but it has 
also some independent meaning as the model for a open spin-glass system. The model 
was introduced in ||, where was shown that random frozen uncorrelated asymmetric 
gaussian bonds (even when the asymmetry is small but generic) totally break finite- 
temperature spin-glass transition. The second important step in this direction was done 
in [ 10 1 . Here multispin interaction spherical spin-glass was considered, and by numerical 
methods was shown that indeed the Langevin dynamics, which is started from the typical 
initial conditions, tends to the paramagnetic state for any nonzero temperature. But if 
the initial conditions are chosen in some special way, then finite-temperature spin-glass 
phase transition to a metastable state is possible. Gugliandolo et. al. [JTUJ proposed 
that this difference between two-spin and multispin interaction models arises due to the 
different structure of the phase spaces when detailed balance holds (i.e., the bonds are 
symmetric). Indeed in the two-spin interaction model the spin-glass phase is marginally 
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stable, hence there is a possibility to destroy the phase transition. In contrast, the 



multispin interaction model has low-lying totally stable states |Tl], [12|, [14 . 

In our model the coupling constants interact with the spins and with a thermal bath at 
the temperature Tj. Generally speaking, can be Tj ^ T a (= T). There are several reasons 
for introducing interactions of this type. A typical example of a two-temperature system 
is a nonequilibrium electron-proton plasma. Due to the large difference between electron 
and proton masses energy transfer between the two components can be neglected in some 
range of times. Protons and electrons then go to equilibrium independently and can have 
different temperatures. In general, a nonequilibrium many- component systems with slow 
energy transfer between the components can, in some range of times, be described by 
introducing several temperatures. Another reason is that: It is well known that in some 
cases and in some sense neural networks can be described by spin-glass models. The spins 
and coupling constants are identified as neurons and synaptic connections. But from a 
physiological point of view the nonlinear interactions with the neurons, and varying the 
synaptic connections with external environment is a very important part of recognition 
and memory. 

At first, spin-glass models with slow varying bonds were considered by Horner |J Q 
Here the model with the detailed balance condition is considered, and the slow motion 
of bonds is used as some method for obtaining correct long-time limits from dynamical 
equations. Other methods for this purpose are also developed by Horner and co-workers 



rj. 



The Langevin equations for the model have the following form: 

Td t Gi = -roi JijCTj + T]i(t), (r]i(t)r)j(t')) = 2r%5(t - t'), (1) 
j 

TjdtJij = ^{ijOiOj - 2vf3jJ ij + r)ij(t), (r]ij(t)Vi'j'(t')) = 2Tj5u>5jf5(t - t'), (2) 

where the bonds (spins) interact with the thermal bath at temperature Tj (T a ), 7 = 



2/N is the standard norm-factor for mean-field models, r is the lagrangian factor for 



1 We thank referee for pointing these references. 
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the spherical constraint X^of = N (here N is the number of spins), Jy 7^ Jji, and uJ?- 
is the potential energy for a coupling constant. 

If the bonds are symmetric a priori, then Eqs. ([!]) and (H) can be viewed as equations 
generated by the mean-field Hamiltonian 



J 2 

ij ij 



It is the Langevin dynamics of the well-known two-spin spherical model where the (sym- 
metric) bonds interact with the thermal bath at the temperature Tj [13]. In model 



considered in this paper asymmetric and symmetric bonds have the similar potential 
energy, and are subjected to the different thermal histories of the same thermostat. 

The analysis of Eqs. ([TJ) and (^) simplifies considerably in the thermodynamic limit 
when N — > 00, where the dynamics of the system can be described by a self-consistent 
equation involving a single spin only. This is achieved by introducing generating func- 
tional for the Langevin dynamics fTS] . Because our interest is in the dynamics of the 



spins only, the integration by the bonds can be taken in the generating functional. 
The resulting mean- field equations read 

{rd t + r)a(t) = fi f dte~ a ^C(t, t)a(t) + rj(t), 
J —00 

(r}{t)r}{t')) = 2r8{t-t')+iiexp(-a\t-t'\)C{t,t'), (3) 

where a = 2v/3j/tj, p, = 2(3(5j/tj, /1 = j3 2 /v/3j. Our purpose here is to get equations 
for the correlation function C(t,t') = (o~(t)o~(t')) and the response function G(t,t') = 
(da(t)/d/3h(t')). Because the initial time — > —00 (it is one of the limits of integration in 
Eq. dpj)) the equilibrium regime is expected 

C(t, t') = C(t - t'), G(t, t') = G(t - t'). 

It should be stressed again that the initial time tends to minus infinity after the thermo- 
dynamic limit. The resulting equations are 

(rd t + r)C(t) = p, I die- a{t -^C{t - t)C(t) + n f°° dte~ a{t+T) C {t + i)G(t) (4) 







(rd t + r)G(t) =fi f dte- a{t - l) C(t - t)G(t) (5) 

Here the difference of the times t — t' is denoted by the same letter t. We also take r = 1 
as fixation of the units. 

Further, we shall investigate the adiabatic limit, where the bonds are frozen with 
respect of the dynamics of spins || Indeed, in a qualitative manner we feel that any 
spin-glass order is possible only if bonds are frozen. This intuitively obvious fact has been 
rederived rigorously recently for spin-glass systems when the detailed balance condition 
holds §, [L3|. Thus the limit 

Tj > t - t' (6) 

is assumed: The bonds are frozen when the spins are observed. Actually, this assump- 
tion is the standard one for investigating the dynamics of usual spin-glass models when 
the bonds are frozen a priori. The different factor of our approach is the noninfinite 
temperature Tj, which is a physical mechanism for introducing a correlation between 
bonds. 

Now if t — t' is large enough, but stills much smaller than tj and all these times are 
much smaller than modulo of the initial time which — > — oo just after N — > oo, we expect 
that the correlation function tends to the Edwards- Anderson order parameter 

C{t - t' h+ oo) = q. 

By the adiabatic condition (|^) in the last integral of Eq. (|j) we can take e~ at ~ 1 
by considering the relevant domain of the integration. In the first integral we have by 
interchanging variable 

jl I" dte- al C p ~ l (t)C(t -i)^—. (7) 
Jo v 

Similar analysis shows that the right-hand term in Eq. (|5|) can be neglected, and we 

have the simple solution 

G(t) = e~ rt (8) 

There are three relevant parameters in our model: v, f3, j3j. In this model we consider the 
phase diagram at the following conditions: v = Tj, n = T/Tj is fixed, and only (3 = 1/T 
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is varying. Thus by 

n i-> (9) 

we go to the model was considered in ||. This fact can be derived directly from (^|): By 
limits (§, the bonds can be considered as uncorrelated, quench, asymmetric gaussian 
variables. 

After this we have the following equation for C(t): 



ip t + r)C(t) = n(3 2 q 2 + (3 2 / dtCit + t)e" ri (10) 

Jo 

This equation holds for t > 0. 

Equation fllUp can be solved by the Laplace transformation. As the general solution 
we have 

= I^V^ + fc+e *v^ + jfe_ e -V^\ (11) 

The constants r should be fixed by the following standard conditions: C(0) = 1 

(the spherical constraint), d t C(t)\ t= o = — 1, A; + = (the condition for monotonic decay 
of C(i)). After some calculations we have for r and Edwards- Anderson parameter q 

n(3 2 q 2 r 



r 2 -f3 2 



(12) 



Finally we have the following equation for q: 



r =\P 2 + 7T^- (13) 
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^Y#(l-!) 2 + l (14) 



1-g 

Besides paramagnetic solution g = there is also the first-order phase transition into 
a spin-glass phase (it is obvious that second-order transition is impossible in this case). 
The transition point is defined as the first temperature when predicts a non-zero 



2 Furthermore, if the bonds in (|l|, ^J) are symmetric a priori then by (^, ^) we go to the usual spherical 
mean- field model [ fL3| . Here n plays the role of "replica number" which in other approaches is introduced 
by the replica trick. 
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(non-paramagnet) solution. The temperature and the jumping of q at this transition 
point can be obtained by the following equations: 

Qc = 3?T2' (15) 

T C = i^, (16) 
where ^ is the positive solution of the equation 

2?/> + 1 = nip(l + tp) 3/2 . (17) 

When n i— > (in this limit we go to the particular case that was studied in Ref. ||) we 
have T c i— > 0, g c i— > 1/3. 

Further, we should discuss the stability properties of the spin-glass solution. If detailed 
balance holds, then static properties of a system is described by the Gibbs distribution 
therefore the stability of solutions can be investigated by analyzing minimums of free 
energy Q. But for this concrete problem free energy is not necessary: Unstability of a 
solution is reflected by unstable behavior of purely dynamical quantities such as positive 
derivation of correlation function, or negative susceptibility. For example, the famous 
(AT) line in usual spin-glass models can be recovered by this purely dynamical consid- 
eration [12, 17]. This property does not connect with any specific character of spin-glass 



systems but is supported by general theorems about the stability of stochastic dynamical 

systems (generally speaking without detailed balance) [Jj, 0, [7j. Because in our case there 

are no such anomalies both for paramagnet and spin-glass solutions we conclude that 

these phases are stable (by chosen initial conditions). The stability of paramagnet for 

any temperature is more or less typical for first-order phase transitions [|ll], |14], [13| . 

Now about an other important problem: If we have two different stable phase, then 

one should be only metastable and another one true stable. In usual models this question 
3 A remark should be given here: Because the mean-field dynamical equations are obtained by N — > oo 
before the initial time — > — oo, some states that arise in dynamics cannot be reflected in the purely static 
investigation by means of free energy. Thus, phase transitions that are predicted by statics and dynamics 
can be different f|, 
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is investigated by free energy: For the true stable phase the free energy should be minimal. 
In principle, a generalized free energy can be introduced also in models without detailed 
balance; the corresponding quantity is based on the so called Kullbak entropy or relative 
entropy |2|. For systems with non-gibbsian stationary distribution this function plays 
nearly the same role as usual free energy but, generally speaking, has a strong dependence 
from initial conditions — in the usual case the detailed balance condition removes this 
dependence. Unfortunately, we have not succeeded in calculation of this function for our 
model. So really we found only that the spin-glass phase occur at least as a metastable 
state. On other hand all stable spin-glass solutions — known by other models (see, for 



example, [|T^, [13], [LB]) — are true stable with respect to paramagnet at least if temperature 
is low enough (of course it is not an argument for true stability in our case but rather 
a hint). One of the scenarios for phase transitions realized in such systems is as follows 
fT3| , |15|| : At some temperature the nontrivial phase (spin-glass in our case) occurs at the 
first time, but true thermodynamical transition occurs at some lower temperature where 
free-energies of the different solutions are the same. May be the scenario for our phase 
transition is the same. 

In the low-temperature limit we have for our solution 

T 2 

1 _ q (18) 

n 

We get this equation by the assumption that 1 — q is small. It should be noted that the 
characteristic time for the correlation function relaxation is t re \ = 1 — q. Thus, critical 
slowing down occur only near T = 0. 

Now about the origin of this spin-glass phase. As we mentioned above behavior of a 
system without detailed balance can be connected with behavior of the phase space of 
the same system when detailed balance holds (i.e., the bonds are symmetric in our case). 
If we assume that such arguments can be used, then at this qualitative level the situation 
is clear: as was shown in [[13] a finite n breaks marginality of the spin-glass solution in the 
two-spin spherical model. 'Replicon' (or dangerous) eigenvalue is ~ n (except exactly the 
critical point). Thus, at least, some spin-glass states are totally stable, and can 'struggle' 
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against the nonhamiltonian influence. 

We consider the spherical spin-glass model with asymmetric partially annealed bonds 
and show that this last property can induce a transition to the spin-glass phase. Several 
very important problems in this field should be considered in future: The problems of 
proving true stability for the new phase, considering partially asymmetric bonds, and of 
course investigating of more realistic models. 
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